Phase transitions in simplified models with long-range 

interactions 



T. M. Rocha Filho, M. A. Amato, B. A. Mello and A. Figueiredo 

Instituto de Fisica and International Center for Condensed Matter Physics 
Universidade de Brasilia, CP: 04455, 70919-970 - Brasilia, Brazil 
(Dated: January 19, 2013) 

Abstract 

We study the origin of phase transitions in some simplified models with long range interactions. 
For the self-gravitating ring model, we were unable to observe a possible new phase transition 
predicted in a recent paper by Nardini and Casetti from an energy landscape analysis. Instead of 
it we observed a sharp, although without any non-analyticity, change from a core-halo to a core 
only configuration in the spatial distribution functions for low energies. By introducing a new 
class of solvable simplified models without any critical points in the potential energy, we show 
that a behaviour similar to the thermodynamics of the ring model is obtained, with a first order 
phase transition from an almost homogeneous high energy phase to a clustered phase, and the 
same core-halo to core configuration transition at lower energies. We discuss the origin of these 
features for the simplified models, and show that the first order phase transition comes from the 
maximization of the entropy of the system as a function of energy and an order parameter, as 
previously discussed by Hahn and Kastner, which seems to be the main mechanism causing phase 
transitions in long-range interacting systems. 

PACS numbers: 05.70.Fh, 05.20.-y, 95.10.Ce 
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I. INTRODUCTION 



The Statistical Mechanics of systems with long-range interactions presents many chal- 
lenging problems, well illustrated by the existence of quasi-stationary non-Gaussian states 
with diverging life-times as the number of particles increases, negative microcanonical heat 
capacity, inequivalence of ensembles and non-ergodicity jl-6]. Examples of systems with 
.ong-range forces cornpr.se self-gravitating systems QB, rron-neutral plasmas flQ, two- 
dimensional flows and some models such as the free electron laser [12J, wave-particle 
interaction in plasmas [l3| and the self-gravitating ring model An interaction po- 

tential is long ranged if it decays for long distances as r~ u , with v < d, and d the spatial 
dimension. Since different statistical ensembles may not be equivalent for long-range forces, 
the choice of the working ensemble must be dictated by the physical constraints on the 
system, and the microcanonical ensemble is the fundamental ensemble in the sense that 
information for other ensembles can be obtained from the knowledge of the microcanonical 
entropy, the converse being not always true [16 |. As a consequence the nature of phase 
transitions may also depend on the physical constraints. It is interesting to note that phase 
transitions occur in one-dimensional classical systems with short-range interactions provided 



some conditions are satisfied 

monly observed Q, Q, Q . 



17H19}. while for long-range interactions they are more com- 



The study of self-gravitating systems with a large number of particles (as in globular 
clusters or in galaxies) is a very difficult task, and the determination of at least some of its 
properties from statistical mechanics is one direction that have been explored during the 
last decades (see 21| and references therein). Nevertheless some problems are still open 
and simplified models have been introduced to clarify those issues such as the determination 



of the distribution function after the process of violent relaxatio n |22H 24|. the dynamics of 
gravothermal collapse [25}, and the nature of phase transitions [26[. It is known that no 
global maximum for the entropy of a self-gravitating system exists unless it is enclosed in 
a container (usually an spherical one with a given large radius), and the 1/r potential is 
somehow bounded from below by a cutoff, by the introduction of a small softening parameter 
e leading to a potential in l/(r + e), or by considering the Pauli exclusion principle for self- 
gravitating fermions (3, I2I]. Examples of simplified models include the self-gravitational 



sheet model 



271 ]. the ring model [14J and the Hamiltonian Mean Field (HMF) model 20], 



all one-dimensional, and two-dimensional gravity jl, 8|. 

In this paper we are interested in discussing the causes of phase transitions in simplified 
models with long-range interactions. Different approaches are reported in the literature to 



study equilibrium phase transitions in those systems [ll, |26|, |28|, |50[. The energy landscape 
method have been used to predict phase transitions in systems with short-range interactions, 
and in some solvable models with long range interactions, based on the topological approach 



to phase transitions 



29 



3jJ . This approach can be briefly summarized as follows. For a 



classical system described by a Hamiltonian of the form 



(1) 



where N is the number of particles and V(q) is the interaction potential, the topological 
approach relies on the study the topology changes of the following manifolds in the config- 
uration space Tat: 

M v = {qE T N \V(q u ...,q N )< Nv}, 



and 



{qeT N \V( qi ,...,q N )=Nv}, 



(2) 



(3) 



with v a fixed value for the potential energy per particle. For most systems, both manifolds 
lead to the same thermodynamical properties [17] . This is not necessarily the case for long- 
range interacting systems (see discussion in section III). The study proceeds using Morse 
theory, where a smooth potential energy is a Morse function (see 29j-[33j for details). For 
short-range confining potential, Franzosi and Pettini presented a proof for a theorem ensuring 
that a phase transition is always associated to a topology transition in the manifold Ai v at 



V = V c = Nv c the critical value of the potential at the phase transition 3414361] . In fact, 



there is a one-to-one relation between topology changes and critical points q of V(q), such 



that V q V\ 



Qi r 1 9=<7 



0, i 



N, provided V is a proper Morse function 



37] . Such critical 



points are also called saddle points in the literature on energy landscapes. 

On the other hand, not all critical points are associated to phase transitions, and more 
recently a new necessary condition, the so-called Kastner-Schnetz-Schreiber (KSS) criterion, 
was formulated stating that only asymptotically flat critical points can lead to a phase tran- 
sition in the thermodynamic limit (see 37| for more details). The microcanonical entropy 
per particle (or degree of freedom) s(e) = S(E)/N, e = E/N, can be non-analytic in the 



thermodynamic limit at a given energy E = E c if: (i) there is a sequence {q^}n=i of 
critical points of iV-particle potential V such that 



lim V(q {N) ) 

N— >oo 



NVr. 



(4) 



where v c = (V)e=eJN is the statistical average value of the potential energy per particle 
at the transition, and (ii) the Hessian matrix of V at the critical points 7f N \ given by 



H 



(AT) 



' d 2 V 
dqidqj 



(5) 



satisfies 



lim 



Det (<>) 



q=q 



1/N 



(N) 



(6) 



Equation flS]) expresses the asymptotic flatness condition for the critical point. This approach 
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was applied with success to the mean field XY and fc-trigonometric models 

Some systems nevertheless are known to display phase transitions not related to a topol- 
ogy change in the manifold Ai v , and particularly in some cases with long-range interacting 
potentials [30[. On the other hand, the KSS criterion was recently applied to the self- 
gravitating ring model, indicating the possibility of a new phase transition hitherto unde- 



tected 



39| . By going beyond previous studies of this model [14j, |l5j, we were unable to 



observe the predicted new phase transition. Instead, we observed a structural transforma- 
tion at lower energies from a core-halo to a core only configuration, such that the halo 
rapidly drops to almost zero. It is reminiscent of the structural change discussed by Aron- 
son and Hansen for an isothermal sphere with a short distance cutoff j^oj (see also In 
order to further discuss these issues, we introduce a new class of solvable simplified models 
with attractive long-range forces with the same qualitative behaviour as the ring model. We 
also discuss for long-range interacting systems admitting a mean-field description a possible 
common origin of phase-transitions. 

The structure of the paper is as follows: in section II we discuss the thermodynamics of 
the ring model analyzing the conjecture by Nardini and Casetti on the possible existence of 
an additional phase transition. In section III we introduce a class of solvable models, and 
discuss with more detail a representative system of this class. We close the paper with a 
discussion of our main results and some concluding remarks in section IV. 



II. THERMODYNAMICS OF THE RING MODEL 



The ring model consists of a system with N particles of mass m constrained to a circular 
ring of radius R, and interacting by their mutual gravitational potential. By a judicious 
choice of units, the Hamiltonian for the model is written as jl^j : 

-i N i N 

where the pair interaction potential: 

Vij = i 1 (8) 

V2y/l - cos(0i -Oj) + e 

with e > the softening parameter, introduced in numerical simulations of self-gravitating 



systems in order to avoid the divergence of the potential for r — > [41]. The coordinate 
9i G [— 7r, 7r) stands for the angle position of the z-th particle on the circle and the 1/N Kac 
prefactor in the second terms in the right hand side of eq. (JTj) is introduced so that the total 
energy of the system is extensive, albeit not additive, and can be achieved by a choice of 
time units. The minimum possible energy per particle is 6q = — l/(2y2e), and in this work 
we rescale the energy unit such that the minimal energy per particle is —1 for all values of 
e. This results in the pair interaction potential: 

Vii = V(6 t - 6j) = — 2 f . (9) 

sj\ - cos^ - Bj) + e 

This model is known to have a phase transition from a homogeneous phase at higher 
energies to an inhomogeneous phase at lower energies. This transition can be continuous or 
discontinuous for higher or lower values of e respectively. For smaller values of e the system 
also has an energy range with a negative heat capacity. These results can be obtained 



by numerical simulations 



14j or a variational method 15| in the microcanonical ensemble. 



Microcanonical Monte Carlo simulations can also be used and will be explored in this work 



together with a variational method. In Reference [39( Nardini and Casetti applied the energy 
landscape method by identifying different kinds of critical (saddle) points of the potential 
(fixing #i = such that the fixed points are isolated): — ir critical points with N n particle 
at 9 = it and iV — N n particles at 9 = 0, and polygonal critical points with N/p particles 
in the p vertices of a regular polygon. They shown that the sequence of polygonal critical 
points satisfy the KSS criterion for the value of v c = V/N corresponding to the transition 



from the homogeneous to inhomogeneous phases, and also that — it critical points satisfy 
the KSS criterion for critical values of the potential per particle given by: 



_ 8 + 2e[6 + 6(5 + 2 e )] 

[(2 + e) 3 /2 + e 3/2]2 > W 

where this expression is already written in units of |eo|. For almost all values of e the 
transition occurs at very low energies (near e = —1), except for e close to 1. Therefore these 
authors conjectured on the existence of a previously unobserved phase at lower energies, 



with the fraction of particles crossing at 6 = ir for v < v c vanishing (see 39] for details). In 
Reference j^J these same authors introduced a simplified model of self-gravitating particles 
derived from approximations on the ring model. They found that particles form a very dense 
cluster, or core, surrounded by particles in a cloud, or halo structure, moving almost freely. 
For very small energies the fraction of particles in the halo is very small, but never vanishes. 



The model in 42j does indeed reproduce quite closely the thermodynamics of the ring model 



43|or Molecular Dynamics (MD) simulations using 



and has many similarities with the Thirring model [45] . In fact, as we argue below, both 
the ring model and the simplified models in section IV show the same qualitative behaviour 
expected to occur for attractive long-range and regularized potentials. 

Figure [1] shows the caloric curves for some values of e obtained using the Entropy Max- 
imization (EM) method described in Reference [15J]. In this method the distribution is 
represented in a suitable numerical grid and is determined iteratively, ensuring convergence 
of the entropy to its maximum for a given energy in the N — > 00 limit, within a mean- field 
description which is exact in this limit. The same curves can be obtained from a microcanon- 
ical Monte Carlo (MC) as described in 
the fourth-order symplectic integrator in j44j , as shown in the right panel of figure HJ From 
our results obtained using MC and MD simulations, and the EM method, as shown in fig- 
ure [H no sign of a phase transition is observed at the predicted critical value v c ~ —.99935. 
Figure [2] shows the entropy and its first and second derivatives as a function of energy 
obtained using the EM method, for an energy interval comprising the phase transition in 
eq. (jTOl) . and no signature of a phase transition is observed. 

A more interesting feature of this model, also shared by the class of simplified models 
in next section, is visible if one considers the first moments of the spatial distribution p{6). 
Since the latter is an even function all odd moments vanish. We consider here the dispersion 



FIG. 1. Left panel: Caloric curves for the potential in ([9]) for some values of e. Right panel: 
Potential energy as a function of total energy (per particle) for e = on an energy interval 

comprising the value v c ~ —.99935 (dotted line) from the EM (diamonds), MC (circles) and MD 
with N = 100 (crosses) u 
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FIG. 2. Entropy and its first and second derivatives with respect to energy around the expected 
phase-transition in eq. [TUJ The first derivative s'(e) is divided by 100 and the second derivative 
s"(e) by 10 6 such that they fit in the same graphic. 

a of the distribution 

a 2 = (6 2 ) = r d6 6 2 p(6), (11) 

J— IT 

and the fourth momentum (kurtosis) M4 of the reduced variable <p = 9/a: 

M 4 = (0 4 ) = ^ d66 4 p(9). (12) 

The kurtosis is very sensitive to minor changes in the distribution function, and any phase 
transition should come about as a discontinuity in M4 as a function of energy, or its derivative 




FIG. 3. Left panel: Standard deviation a as a function of the energy e for a few values of e. Right 
panel: Fourth moment M4 (kurtosis). 



(and a also). Figure [3] shows a and M4 as a function of energy obtained using the EM 
method. These results were tested for numerical convergence with a maximum relative error 
in the entropy of 10~ 8 at each energy value. For small values of e a very strong maximum for 
M 4 is obtained alongside a fast variation of the derivative of a. This corresponds to a change 
of regime in the distribution function, from a core-halo structure to a core only regime, in 
the sense that in the latter case the halo is so extremely diffuse that it can be considered 
as nonexistent for any practical purposes. The distribution in the halo cannot vanish for 
any coordinate if the mean-field self-consistent description is valid, as it is the solution of 
a self-consistent equation of the form p{9) oc exp(— /3v (0)), and therefore this would imply 
a diverging positive potential. The maximum of M 4 and the corresponding value of the 
potential per particle v as a function of e obtained using the EM method are shown in 
figure H] for values from e = 10~ 26 up to e = 10. From our results it seems that a genuine 
phase transition only occurs for e — > 0, corresponding to a collapse of the system into a single 
point for T = 0. The corresponding values of the predicted phase transition in eq. ([TO]) are 
also shown in the right panel of fig. S] and exhibits a similar qualitative behaviour as the 
curve corresponding to the actual regime change. Even though the presence of a critical 
point satisfying the KSS criterion may somehow induce this change of regime, we will show 
below that simplified models with no related critical points can exhibit the same type of 
behaviour. 
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FIG. 4. Left panel: Maximum vale of M4 as a function of e. Right panel: Value of the potential 
per particle at the Kurtosis maximum (diamonds) and expected phase transition from eq. (jlOp 
(continuous line) function of e. 



III. A CLASS OF SIMPLIFIED MODELS 



In this section we discuss a class of solvable models with long range interactions, exhibiting 
a behaviour in many aspects similar to the ring model, with a first order phase transition 
from a quasi-homogeneous to a strongly inhomogeneous states, despite the fact that its 
potential energy has no critical (saddle) points. Let us consider a pair interaction potential 
of the form: 

Vij = -uir^uijj), (13) 

with u(r) a smooth positive function of the distance to the origin of the coordinate system 
such that Vij is attractive, and the Hamiltonian given by eq. ([7]). The number of dimensions 
is restricted here to one, but similar solvable models can be obtained for any dimensions 
following the same procedure. The class of models corresponding to eq. ( 1T3|) comprises both 
the Thirring (with a Kac 1/N regularization factor) and HMF (at equilibrium) models as 



particular cases 



45 



46| . This whole class of models is solvable in the limit iV — > 00 (see 



below). The mean field description is exact in this limit and the iV-particle equilibrium 
distribution function fully factorizes as j^ : 

N 



f N (pi,r 1 ,...,p N ,r N ) = Y[fi{pi,ri) + 0(1/N), 



(14) 



where fi is the one particle (equilibrium) distribution function. The equilibrium properties 



can be computed using large deviation techniques 



48 



49], but here we chose to use the 



approach in 50[ which is more intuitive for the present purposes. The equilibrium properties 
of the system in the microcanonical ensemble can be determined by maximizing either the 
Boltzmann or Gibbs entropy 5jj, respectively: 



S B (E) = log 



J 5{H{p,q)-E) dT 



(15) 



S G (E) = -J f N lnf N dT, (16) 

where H (p, q) is the iV-particle Hamiltonian with canonical variables collectively denoted by 
p, q, dT a volume element in phase space and the Boltzmann constant is set to unity. The 
Gibbs entropy is more appropriate in the present context as it will become clear below. 
For a fully factorized distribution function, the Gibbs entropy per particle is written as 



s e 



dpdr fi(p,r) ln/i(p,r). 



(17) 



The equilibrium distribution is obtained by maximizing s(e) with the normalization con- 
straint: 

J dpdrfx{p,r) = 1, (18) 

and the total energy constraint: 

I dp dr ^—fi(p,r) — - f dpdr dp' dr' u(r)u(r')fi(p,r)fi(p',r') = E, (19) 

with m the particle mass. Local maxima and saddles of the entropy correspond to meta- 
stable and unstable states, respectively. We therefore obtain: 



SS = -\nf 1 -l-\-(3e(p,r), 



(20) 



with A and Lagrange multipliers corresponding to constrains ( Tl8|) and (fl9|) respectively, 
the mean- field energy per particle: 

,2 



P 

e (P; r ) = 7T~ + U ( r ) 
2m 



(21) 



and the mean-field potential: 

U(r) = J dr'V(r,r') = -u(r) J dp' dr 1 'u{r')f x {p' >') = -cra(r), (22) 



with 



a 



J dp' dr' u{r')fi{p\r'). 



(23) 
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Equation (|20|) then implies that the one particle distribution function is the solution of the 
self-consistent equation: 



fi(p,r) = Cexp 



B ( a u(r) — ^~ 
' \ K J 2m, 



with C = exp(— A— 1) the normalization constant. Plugging eq. 
over the moments we finally obtain: 



into 



a /I3= f dru{r)e au{r) x 



dr e au{r) 



(24) 

and integrating 
(25) 



with a = Pa. Equation f )25|) is a transcendental equation for a (or a), and the equilibrium 
distribution function is then obtained from eq. ( 124")) . The stability of this solution must also 



be determined, and may depend on the ensemble considered if they are not equivalent 



48 



50], i. e. we must impose that f\ maximizes the entropy or minimizes the free energy, for 



the microcanonical and canonical ensemble, respectively. Different models are obtained by 
specifying the physical space range and the function u[r). 
Let us consider the function in one- dimension: 



u{r) 



(26) 



y/r + e 

with r = |sc|, e > a softening parameter and x the spatial coordinate. The system is 
enclosed in a one-dimensional box defined by L\ < x < L 2 and m — 1. The numerator in 
the right-hand side of eq. ( 1261) is chosen such that the minimal energy per particle is —1. In 
order to straightforwardly apply Morse theory the potential must be a Morse function, and 
therefore it must be a smooth (analytical) function in all configuration space. This is ensured 
by taking Li = and L 2 = L > 0. Even though the resulting model is quite artificial, it 
still allows to discuss the main results of this section. The algebraic equation (|25|) is here: 



a//3 = Q(a), 



(27) 



with 



4eEi (-Aa) - 4eEi (-^2 a) + 2v / 2~iv 'L + e e Aa - 2eV2e^ a 

Ei (Aa) - Ei (-v^ct)]}^ 



Q(a) 

x {(L + e)e m - ee^ (y/2a + 1 
where Ei is the exponential integral function and A 



2ea 2 



(28) 



\J2e/\jL + e. The energy per particle 
is obtained from eqs. (1211 and (|22|) and the solution of eq. (1271) : 



1 



a 
~2 



(29) 
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In this expression the potential energy is not the average value of the mean-field potential, 
but half this value. 



Without loss of generality we can choose L = 1. The graph of Q(a) for e = 10 6 is given 
in fig. |5j Equation ( |27|) has three different solutions for 10.6 < f3 < 1386.2, or equivalently 
0.00072 < T < 0.094 and -0.738 < e < 0.00035. Outside this interval the algebraic 
equation (|27|) has a single solution. These different branches can be identified in the caloric 
curve shown in fig. |5j where there are either one or three different energy values for a given 
temperature, corresponding to the possible solutions of eq. (1271) . From figure [6] we see that 
a first order phase transition occurs at e ~ 0.00128 from an almost homogeneous state with 
positive heat capacity to a strongly inhomogeneous state with negative heat capacity. The 
heat capacity changes sign at a lower energy. A measure of the degree of inhomogeneity in 
the high energy phase is given by Q(a = 0), which at its turn is a monotonously increasing 
function of e. In the canonical ensemble the phase transition occurs with an energy jump jlj]. 
This behaviour is qualitatively the same for all small values of e. For e > 1.05 x 10~ 4 there 
is no phase transition. 

The phase transition do not correspond to any stationary point in the potential of the 
system, as it has no stationary points. The only topology change of Ai v occurs at the 
potential minimum corresponding to v = —1, as for v < — 1 the set set hA v is empty. 
This topology change is not associated to any finite temperature phase transition. On the 
other hand, the manifold E„ changes topology at v — — 1 and v — e/(L + e) = v max , the 
latter corresponding to the maximum of the potential, the set E„ being non-empty only for 
— 1 < v < v max . It is noteworthy that the phase transition in this model, as for the ring and 
HMF models, are a consequence of the potential being bounded from above. Therefore one 
might expect that the phase transition is triggered by the topology change of E„ at v — v max . 
Nevertheless the critical value v c of the potential at the transition is always smaller that v max 
as shown in fig. [71 Therefore, no topology change of either M. v or can be associated to 
the phase transition observed in this model. A similar behaviour was also obtained for the 



spherical model [52l. 153]. 



The origin of the phase transition can be explained by considering the entropy as a 
function of both the energy e and the mean-field variable a, obtained from eq. f fl6|) and the 
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-0.4 

e 



FIG. 5. Left panel: Q(a) for e = 10~ 6 , with Q(a = 0) « 0.028. Right panel: Caloric curve. The 
temperature jump is visible on the insert. The dotted line represents the meta-stable and unstable 
branches. 

distribution function in eq. $M§, after eliminating /3 from eq. fl29|) . as: 

s(e,a) = 2/3e- log C- -. (30) 



Note that the normalization factor C in eq. ( |24|) is also a function of e and a. The function 
s(e, a) is shown on the left panel of figure |8]in a region of the (e, a) containing the first-order 
phase transition, where it is clear that it is not always a concave function for all values of 
its arguments. The equilibrium entropy, as obtained above, is then given by: 

s(e) = max Q s(e, a). (31) 



Even though s(e, a) is an analytical function of both its arguments, the maximization in ( )3T|) 
with respect to a yields a function with a non-analyticity at the critical energy, as shown 
in fig. (E]). Hahn and Kastner discussed t 



has a role similar to the variable a 



54 



;his scenario for a system where the magnetization 
551 ] . It is straightforward to see that provided u(r) is 
a positive monotonously increasing function of r, the qualitative behaviour of the resulting 
models is the same as the representative example just discussed above. 

IV. DISCUSSION AND CONCLUDING REMARKS 



We have shown that the conjectured new phase transition for the ring model, as discussed 
in Ref. 39] is not observed in our results, obtained from MC, MD and EM methods. The 
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FIG. 6. Equilibrium entropy S in the first order transition region. The dotted line shows the 
meta-stable and unstable branches. 



FIG. 7. Potential per particle v c at the phase transition for some values of e obtained from the EM 
method, and the maximum value of the potential per particle v max . 



same model nevertheless displays a marked change of regime from a core-halo to a core only 
structure. For the core structure, the probability of a particle being in the halo is so small 
that it can be considered as vanishing for any practical purposes. In the previous section we 
introduced a class of simplified solvable models, and studied in detail a representative case 
as defined by the choice of potential energy in eqs. f JT3|) and ff26|) . This particular model 
has properties very similar to those of the ring and Thirring models, and to the simplified 



model in Ref . 42] . It has a first order phase transition from an almost homogeneous phase 



to a clustered phase, and a core-halo to core structural transition. The important point 
for this class of models, with a monotonous function u(r), is that it is possible to chose a 
potential with no stationary points, such that a first order phase transition is obtained from 
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FIG. 8. Left panel: Entropy as a function of the energy e and mean-field variable a. on the region 
of the first order phase transition. Right panel: Entropy as a function of a for different energy 
values. For some energy values a local maximum exists corresponding to meta-stable states near 
the phase transition 



the maximization of the entropy per particle s(e, a) with respect to the variable a given 
in eq. ( 123]) . Even with stationary points in the potential, the entropy is obtained by the 
maximization in eq. (|3T|) of a smooth non strictly concave function. This origin for phase 



30, 



54 



55] 



47] is such 



transitions in long-range interacting systems was discussed by Hahn and Kastner 
In fact, it seems that the mean- field nature of long range interacting systems 
that this mechanism is ubiquitous in these systems, regardless of the existence of a topology 
change in phase space. As a perspective, it would be interesting to investigate if a similar 
behaviour is also observed for two and three-dimensional self-gravitating systems, as they 
also display a core-halo equilibrium structure js]. 
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FIG. 9. Left panel: Mono-log plot of the dispersion a as a function of energy for e = 10" 6 . Right 
panel: Zoom over the first-order transition region showing the solution for the stable (continuous 
lines), meta-stable and unstable branches (dotted line). 




FIG. 10. Left panel: Fourth momentum M4 of the reduced variable r = r/a. Right panel: M4 in 
the first order transition with the meta-stable and unstable branches given by the dotted line. 
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